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1^-^ ' Abstract 

■ In the present work we establish a simple relation between the Dirac 

| equation with a scalar and an electromagnetic potentials in a two-dimensional 

Q j. case and a pair of decoupled Vekua equations. In general these Vekua 

I ■ equations are bicomplex. However we show that the whole theory of 

| pseudoanalytic functions without modifications can be applied to these 

equations under a certain not restrictive condition. As an example we 
C ' formulate the similarity principle which is the central reason why a pseu- 

. . , doanalytic function and as a consequence a spinor field depending on two 

space variables share many of the properties of analytic functions. 

One of the surprising consequences of the established relation with 
pseudoanalytic functions consists in the following result. Consider the 
Dirac equation with a scalar potential depending on one variable with fixed 
energy and mass. In general this equation cannot be solved explicitly even 
if one looks for wave functions of one variable. Nevertheless for such Dirac 
equation we obtain an algorithmically simple procedure for constructing 
in explicit form a complete system of exact solutions (depending on two 
variables). These solutions generalize the system of powers l,z,z 2 , . . . in 
complex analysis and are called formal powers. With their aid any regular 
solution of the Dirac equation can be represented by its Taylor series in 
formal powers. 

1 Introduction 

The Dirac equation with a fixed energy and the Vekua equation describing pseu- 
doanalytic (=generalized analytic) functions both are first order elliptic systems, 
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and it would be quite natural to expect a deep interrelation between their the- 
ories especially in the case when all potentials and wave functions in the Dirac 
equation depend on two space variables only. Nevertheless there is not much 
work done in this direction 1 due to the fact that any of the traditional matrix 
representations of the Dirac operator does not allow to visualize a relation be- 
tween the Dirac equation in the two-dimensional case and the Vekua equation. 
The Dirac equation is a system of four complex equations which does not de- 
couple in a two-dimensional situation but decouples in the one-dimensional case 
only. 

In the present work we establish a simple relation between the Dirac equa- 
tion with a scalar and an electromagnetic potentials in a two-dimensional case 
from one side and a pair of decoupled Vekua equations from the other. As a 
first step we use the matrix transformation proposed in [S] (see also ^U] an d 
|14| ) which allows us to rewrite the Dirac equation in a covariant form as a 
biquaternionic equation. This is not our aim to discuss here the advantages of 
our biquaternionic reformulation of the Dirac equation compared with other its 
representations (the interested reader can find some of the arguments in |lt)|h 
We point out only that our transformation is C-linear as well as the resulting 
Dirac operator, which is not the case for a better known biquaternionic refor- 
mulation of the Dirac operator introduced by C. Lanczos in JS] (see and JU| 
for more references). Moreover, in the time-dependent case, with a vanishing 
electromagnetic potential our Dirac operator is real quaternionic. 

Here we exploit another attractive facet of our biquaternionic Dirac equation. 
In the two-dimensional case it decouples into two separate Vekua equations. In 
general these Vekua equations are bicomplex. However we show that the whole 
theory of pseudoanalytic functions without modifications can be applied to these 
equations under a certain not restrictive condition. As an example we formulate 
the similarity principle which is the central reason why a pseudoanalytic function 
and as a consequence a spinor field depending on two space variables share 
many of the properties of analytic functions; e.g., they are either identically 
zero or have isolated zeros. In this way more results of the theory developed in 
[TT?| and in posterior works (see, e.g., and can be applied to the two- 
dimensional Dirac equation with a scalar and an electromagnetic potentials. 
Nevertheless in the present work we concentrate on another non-trivial and 
surprising consequence of the established relation with pseudoanalytic functions. 
Consider the Dirac equation with a scalar potential depending on one variable 
with fixed energy and mass. In general this equation cannot be solved explicitly 
even if one looks for wave functions of one variable. None the less the result of 
this work is an algorithmically simple procedure for obtaining in explicit form 
a complete system of exact solutions depending on two variables for such Dirac 
equation. This system of solutions is a generalization of the system of powers 
1, z, z 2 , . . . in complex analysis and as such they are not appropriate for studying 
the Dirac equation on the whole plane. However the very fact that it is always 

1 We refer to the work [T| where the theory of pseudoanalytic functions was used in a way 
completely different from ours for studying the two-dimensional Dirac equation with a scalar 
or a pseudoscalar potential. 
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possible to obtain explicitly a complete system of exact solutions of the Dirac 
equation with scalar potential of one variable as well as the hope to be able to 
obtain explicitly not only the generalizations of positive powers but also those of 
the negative ones makes in our opinion this approach attractive and promising. 
The system of exact solutions for the Dirac equation with a one-dimensional 
scalar potential is obtained due to the proposed reduction of the Dirac equation 
to Vekua equations and due to L. Bers' theory of Taylor series in formal powers. 

In Section 2 we introduce notations. In Section 3 we give the biquaternionic 
reformulation of the Dirac equation. Let us emphasize that our biquaternionic 
Dirac equation is completely equivalent to the "traditional" Dirac equation writ- 
ten in 7-matrices, we have a simple matrix transformation giving us a relation 
between their solutions. In Section 4 we show that in a two-dimensional situation 
the Dirac equation with a scalar and an electromagnetic potentials decouples 
into a pair of bicomplex Vekua equations. We establish that if one of the co- 
efficients in such Vekua equation has not zeros and does not turn into a zero 
divisor at any point of the domain of interest, the solutions will not be zero 
divisors either, and the whole theory of generalized analytic functions without 
modifications is applicable to the bicomplex Vekua equation. 

In Section 5 we adapt some definitions and results from L. Bers' theory to 
bicomplex pseudoanalytic functions. Section 6 is dedicated to a special class 
of Vekua equations which have been studied recently (see m and [T2]) due to 
their close relation to stationary Schrodinger equations. In Section 7 we show 
that the Dirac equation with a scalar potential depending on one space variable 
can be represented as a Vekua equation from the special class mentioned above. 
Here we should notice that the case of the scalar potential is only an example. 
The same is true, for example, for the electric potential. To the Vekua equation 
we apply L. Bers' procedure for constructing corresponding formal powers which 
as was mentioned above are exact solutions of the Vekua equation and generalize 
the system of analytic functions 1, 2, z 2 , . . .. With their aid any regular solution 
of the Vekua equation can be represented by its Taylor series in formal powers. 

2 Preliminaries 

We denote by H(C) the algebra of complex quaternions (= biquatcrnions). The 
elements of H(C) have the form q — J2k=o Ik^k where {qk} C C, cq is the unit 
and {efc|fc = 1, 2, 3} are the standard quaternionic imaginary units. 

We denote the imaginary unit in C by i as usual. By definition i commutes 
with efe, k — 0,3. We will use also the vector representation of q 6 H(C) : 
q = Sc(q) + Vec(q), where Sc(q) = q a and Vec(q) = ~q = ^ fc=1 g^e*,. The 
quaternionic conjugation is defined as follows q = q — ~q . 

By M p we denote the operator of multiplication by p from the right hand 
side 

M p q = q-p. 

More information on complex quaternions the interested reader can find, e.g., 
in OH or [Tlj. 
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Let q be a complex quaternion valued differentiable function of x = x%, 2:3). 
Denote 



This operator is called sometimes the Moisil-Theodorescu operator or the Dirac 
operator but the true is that it was introduced already by W.R. Hamilton himself 
and studied in a great number of works (see, e.g., (SJ, 0, [El, (TTij . 52)- 

3 Quaternionic reformulation of the Dirac equa- 
tion 

Consider the Dirac operator with scalar and electromagnetic potentials 



where 7-,-, j — 0, 1, 2, 3 are usual 7-matrices (see, e.g., @], m G IR, p e i, Ak 
and p sc are real valued functions. 

In |H] a simple matrix transformation was obtained which allows us to rewrite 
the classical Dirac equation in quaternionic terms. 

Let us introduce an auxiliary notation / := f(t, x\, x%, — X3). The domain 
G is assumed to be obtained from the domain G C I 4 by the reflection X3 — > 
— X3. The transformation announced above we denote as A and define it in the 
following way. A function $ : G C M 4 — > C 4 is transformed into a function 
F:GcI 4 ^ H(C) by the rule 



Let us present the introduced transformations in a more explicit matrix form 
which relates the components of a C 4 -valued function $ with the components 
of an H(C)-valued function F: 






The inverse transformation A 1 is defined as follows 



$ = A-^F] = (-iFx - F 2 , -F - iF 3 ,F - i%,iF x - F 2 f '. 
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Denote 

R = D — d t M ei + a + A/-*(p^i-*(p-+>n)e 2 ) 
where a = i(A\e\ + A 2 e 2 — ^363). The following equality holds [TU] 

i? = ^7 1 7 2 7 3 D^- 1 . 
That is, a C 4 -valued function $ is a solution of the equation 

D$ = in G 

iff the complex quaternionic function F = AQ is a solution of the quaternionic 
equation 

RF = in G. 

Note that in the absence of the electromagnetic potential the operator R be- 
comes real quaternionic which is an important property (see |13|). 

In what follows we assume that potentials are time- independent and consider 
solutions with fixed energy: x) = $ w (x)e lwt . The equation for $ w has the 
form 

=0 in G (1) 



where G is a domain in M 3 , 

3 / 3 \ 

EL = iw7o + ^ "fkdk +i[m + p e aa + ^ Ak ^ k +Psc\ ■ 

k=l V fc=l / 

We have 

R u = ^7i7273D^^ _1 , 

where 

R u = D + el + M h 

with b = — i((p e l + w)ej — i(p« c + m)e 2 ). Thus, equation QJ turns into the 
complex quaternionic equation 

R^q = (2) 
where q is a complex quaternion valued function. 
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4 The Dirac equation in a two-dimensional case 
as a bicomplex Vekua equation 

Let us introduce the following notation. For any complex quaternion q we denote 
by Qi and Q2 its bicomplex components: 

<2i = go + <?3e3 and Q 2 = q2 - qie 3 - 

Then q can be represented as follows q — Qi + Q 2 e 2 - For the operator D we 
have D = Di + D2&2 with D\ — e^d^ and D2 = 82 — d\e 3 . Notice that b = Be2 
with B = —(psc + m) + i(p e i + w)e3, a = A\ + A2&2 with A\ — a^e^ and 
A2 = a 2 - aie 3 . 

We obtain that equation J2J is equivalent to the system 

D1Q1- D2Q2+ A1Q1- A2Q 2 -BQ 2 = 0, (3) 

D2Q 1 + D 1 Q2+A 2 Q 1 +A 1 Q2 + BQ 1 = 0, (4) 

where Q\ and Q2 are bicomplex components of q. We stress that the system 
©, (@J is equivalent to the Dirac equation in 7-matrices Q). 

Let us suppose all fields in our model to be independent of £3, and A\ = 
0363 = 0. Then the system ifflk Q decouples, and we obtain two separate 
bicomplex equations 

D2Q2 = -A2Q2 - BQ 2 

and 

D 2 Qi = -A 2 Qi -BQ V 

Denote d = D 2 , a = —A 2 , b = —B } w = Q 2 , W = Qi, z — x + yk, where 
x = X2, y = xi and for convenience we denote k = e^. Then we reduce the 
Dirac equation with electromagnetic and scalar potentials independent of X3 to 
a pair of Vekua-type equations 

dw = aw + bw (5) 

and 

OW = aW + bW. (6) 

The difference between the bicomplex equations (JjJ , (jHJ and the usual complex 
Vekua equations is revealed if only w or W can take values equal to bicomplex 
zero divisors (otherwise equations JSJ) , iJBJ can be analyzed following Bers- Vekua 
theory [3], EI)- Let us study this possibility with the aid of the following pair 
of projection operators 

P+ = -(l+«k) and P _ = -(l-ik). 

The set of bicomplex zero divisors, that is of nonzero elements q = qo + qik, 
{go, qi} C C such that 

m = (go + gik) (g - gik) = (7) 

we denote by 6. 
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Lemma 1 Let q be a bicomplex number of the form q = qo + qik., {qg, qi} C C. 
Ifqe & then q = 2P+q Q or q = 2P-q . 



Proof. From (7J| it follows that q$ + q\ = which gives us that q\ = ±iqo. 
That is q = <7q(1 + *k) or q — q (l — ik). ■ 

For other results on bicomplex numbers we refer to |16| . 

Let n denote a bounded, simply connected domain in the plane of the vari- 
able z. 

Theorem 2 Let b(z) & U {0}, Vz G fl and w, W be solutions of £5J) and 0) 

respectively. Then w{z) ^ 6 and W(z) ^ &, \/z G Q. 

Proof. Assume that w(z) S & for some z 6 fl. For definiteness let w(z) — 
2P + wq(z). Then from J^J we have 



Applying to this equality wc find that P~b = which is a contradiction. ■ 
Thus, if the coefficient b does not have zeros and does not turn into a zero 
divisor at any point of the domain of interest, the solutions of JS| and © will 
not be zero divisors either and the whole theory of pseudoanalytic functions 
is applicable without changes to the bicomplex equations |JSJ and As an 
example let us formulate one of the main results of the theory, the similarity 
principle which is the basic tool for studying the distribution of zeros and of 
singularities of pseudoanalytic functions as well as boundary value problems 



Theorem 3 Let w be a regular solution of Jj5J) in a domain and let b(z) ^ 
& U {0}, Vz G £L Then the bicomplex function $ = w ■ e h , where 



The proof of this theorem is completely analogous to that given in |T^]. It 
would be interesting to extend this result to the case of b being a zero divisor 
in the whole domain Q or in some points. 

This theorem opens the way to generalize many classical results from theory 
of analytic functions to the case of solutions of equations (JSJl and 10 by analogy 
with ^]. Nevertheless in the present work we prefer to explore another possi- 
bility. Namely, we show how the application of Bers' theory of pseudoanalytic 
functions allows us to obtain explicitly a complete system of solutions of the 
Dirac equation with a scalar potential depending on one variable. 



dP + w a = aP+wo + bP~w . 





is a solution of the equation 9$ = in Q. 
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5 Some definitions and results from Bers' the- 
ory for bicomplex pseudoanalytic functions 

5.1 Generating pair, derivative and antiderivative 

Following [5] we introduce the notion of a bicomplex generating pair. 

Definition 4 A pair of bicomplex functions F — Fq + i*ik and G — Gq + Gik, 

possessing in 17 partial derivatives with respect to the real variables x and y is 
said to be a generating pair if it satisfies the inequality 

Vec(FG) ^ infl. 

The following expressions are called characteristic coefficients of the pair (F, G) 

FG- — F-G FG- — F-G 



a(F,G) -= =pr, { F,G) 



FG-FG y ' ' FG-FG 



FG, — F Z G „ FG Z — F Z G 



A {F,G) ~ T= ^7T' B (.F,G) 



FG-FG v ' ' FG-FG 

where the subindex z or z means the application of d or d respectively. 

Every bicomplex function W defined in a subdomain of f2 admits the unique 
representation W = <j)F + tpG where the functions <f> and ip are complex valued. 

The (F, G)-derivative W = d(F '^ z — of a function W exists and has the form 

W = cj> z F + ^ z G = W z -A {FjG) W-B {FtG) W (8) 

if and only if 

<fcF + tp-G = 0. (9) 
This last equation can be rewritten in the following form 

which we call the bicomplex Vekua equation. Solutions of this equation are 
called (F, G) -pseudoanalytic functions. 

Remark 5 The functions F and G are (F, G)-pseudoanalytic, and F = G = 0. 

Definition 6 Let (F,G) and (F±,Gi) - be two generating pairs in 57. (Fi,Gi) 
is called successor of(F,G) and (F,G) is called predecessor o/(iq,Gi) if 

a (fi,Gi) = a (F,G) an d k(F 1: Gi) = —B(F,G)- 

The importance of this definition becomes obvious from the following state- 
ment. 
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Theorem 7 Let W be an {F 1 G)-pseudoanalytic function and let (Fi,G\) be a 
successor of (F, G). Then W is an {F\,Gi)-pseudoanalytic function. 

Definition8 Let(F,G) be a generating pair. Its adjoint generating pair (F, G)* 
(F*, G*) is defined by the formulas 

i^* ^F s-<* 2G 

F = = — =— , G = 



FG-FG FG-FG 
The (F, G)-integral is defined as follows 

J Wd {FiG) z = ^F(z 1 )Sc J G*Wdz + G( Zl )Sc J F*Wdz^j 

where T is a rectifiable curve leading from zq to z\. 

If W — <pF + ipG is an (F, G)-pseudoanalytic function where <f> and ip are 
complex valued functions then 



/' 



Wd (F , G) z = W(z) - <t>(z )F(z) - 1>{zo)G{z), (10) 



and as F =G = 0, this integral is path-independent and represents the (F, G)- 
antiderivative of W. 

5.2 Generating sequences and Taylor series in formal pow- 
ers 

Definition 9 A sequence of generating pairs {(F mi G m )}, m = 0, ±1, ±2, . . . , 
is called a generating sequence if (F m+ i, G m+ i) is a successor of (F m ,G m ). If 
(F ,G ) = (F,G), we say that (F,G) is embedded in {(F mi G m )}. 

Theorem 10 Let (F, G) be a generating pair in 0. Let 17 1 be a bounded do- 
main, fli C ft. Then (F,G) can be embedded in a generating sequence in Q\. 

Definition 11 A generating sequence {(F mi G m )} is said to have period /i > 
if (Fm+fj,, G m+tl ) is equivalent to (F m ,G m ) that is their characteristic coeffi- 
cients coincide. 

Let W be an (F, G)-pseudoanalytic function. Using a generating sequence 
in which (F, G) is embedded we can define the higher derivatives of W by the 
recursion formula 

W^=W; W ^] = d (^f Wlm \ m =l,2,.... 

dz 
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Definition 12 The formal power Z$(a, zq; z) with center at zq 6 fl, coefficient 
a and exponent is defined as the linear combination of the generators F m , G rn 
with complex constant coefficients X, fj, chosen so that XF m (zo) + fj,G m {zo) = a. 
The formal powers with exponents n — 1,2,... are defined by the recursion 
formula 

Z% +1 \a,Z ;z) = (n + 1) f Z^ +1 (a,z ; C)^ m , Gm )C- (H) 

This definition implies the following properties. 

1. Zm\a, zo; z) is an (F m , G m )-pseudoanalytic function of z. 

2. If a' and a" are complex constants, then 

+ ka", z ; z) = a'4"'(l, z ; z) + a"Z^(k, z ; z). 

3. The formal powers satisfy the differential relations 

j- z = nZ m+i ( a > z o; z). 

4. The asymptotic formulas 

Zffl (a, z \ z) ~ a(z - z ) n , z -> z 

hold. 

Assume now that 

oo 

W(z) = £zW(o,2!b;«) (12) 

Tl = 

where the absence of the subindex m means that all the formal powers corre- 
spond to the same generating pair (F,G), and the series converges uniformly 
in some neighborhood of Zq. It can be shown that the uniform limit of pseu- 
doanalytic functions is pseudoanalytic, and that a uniformly convergent series 
of (F, G)-pseudoanalytic functions can be (F, G)-differentiated term by term. 
Hence the function W in (|12l) is (F, G)-pseudoanalytic and its rth derivative 
admits the expansion 

oc 

wW(z) = n(n - 1) • • • (n - r + \)Z^- r \a n , z ; z). 

n=r 

From this the Taylor formulas for the coefficients are obtained 

a. . (13) 

71! 
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Definition 13 Let W(z) be a given (F,G)-pseudoanalytic function defined for 
small values of \z — zq\. The series 

oo 

J2z {n \a,zo;z) (14) 

n=0 

with the coefficients given by flty) is called the Taylor series ofW at zq, formed 
with formal powers. 

The Taylor series always represents the function asymptotically: 

N 

W(z)-Y,Z {n) (a,z ;z)^o{\z-z \ N+i y z - z , (15) 

n=0 

for all N. This implies (since a pseudoanalytic function can not have a zero 
of arbitrarily high order without vanishing identically) that the sequence of 
derivatives { W'"'^)} determines the function W uniquely. 

If the series (|14|) converges uniformly in a neighborhood of zq, it converges 
to the function W. 

Theorem 14 The formal Taylor expansion \14\l of a pseudoanalytic function 
in formal powers defined by a periodic generating sequence converges in some 
neighborhood of the center. 

6 Special class of Vekua equations 

The following important class of Vekua equations was considered in \T2\ . Let 
fo be a complex valued (with respect to i), twice differentiable nonvanishing 
function defined on Q. Consider the equation 

dW = ^fe\Y in n. (16) 
Jo 

Denote v x = A/ // . 

Theorem 15 fWf IfW = Wi + W 2 k is a solution of fWj) then W x = ScW is 
a solution of the stationary Schrddinger equation 

-AWi + P1W1 = in SI (17) 
and W2 = Vec W is a solution of the associated Schrddinger equation 

-AW 2 + v 2 W 2 = inQ. (18) 
where v 2 = 2(df ■ 9/o)//o - v i- 

Moreover, in |12j a simple formula was obtained which allows us for any given 
solution W\ of (nil to construct such a solution W 2 of {THl that W — W\ + W 2 k 
will be a solution of (|16|l generalizing in this way the well known procedure for 
constructing conjugate harmonic functions in complex analysis. 
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7 Dirac equation with a scalar potential 



Let us show that the Dirac equation with a scalar potential depending on one 
real variable reduces to a bicomplex Vekua equation of the form <|16|) • 

Let p sc — p(x) and p e i = 0, A/. = 0, k = 1, 2, 3. Then according to Section 
0| the Dirac equation is equivalent to the pair of bicomplex Vekua equations 

dw = bw (19) 

and 

dW = bW (20) 

with b = p(x) + m — iwk. 

Let /o = e p ( x )+ mx + iu) y ^ w here P is an antiderivative of p. Then we have 

b = df a /f Q . 

Note that due to theorem ll5l if the bicomplex function W is a solution of Q20|l 
then the complex function W\ = Sc W is a solution of the stationary Schrodinger 
equation (|17|l where 

Vx(x) =p'(x) + (p(x) +mf - uj 2 , (21) 

and the function Wi — Vec W is a solution of equation (|18fl where 

!/ 2 (x) = -p'(x) + (p(x) + mf - uj 2 . (22) 

Let us notice that both Schrodinger equations 1171) and (|18fl in this case admit 
separation of variables. Nevertheless this does not imply they can be solved 
explicitly. In general this is not the case. However we will show how using our 
approach and Bers' theory for both of them one can construct in explicit form 
a locally complete system of exact solutions. 

Consider equation H20|) . It is easy to see that the pair of functions 

F = f and G=^ (23) 

Jo 

represents a generating pair for H2()(l . Note that F — e a and G = e _CT k, where 
a = a(x) + (3{y) and a(x) — P(x) + mx, (3(y) — iuoy. For a generating pair of 
such special kind it is easy to construct a successor 0. Let r = — a(x) + (i{y)- 
Then the pair Fi = e T and G\ — e~ T k is a successor of (F,G). Moreover, 
(F, G) is a successor of (F±, G\). Thus, for (F, G) we obtain a complete periodic 
generating sequence of a period 2 in explicit form (for explicitly constructed, 
in general non-periodic generating sequences in a far more general situation we 
refer to [E]). 

The fact that we have a generating sequence in explicit form implies that we 
are able to construct the corresponding formal powers of any order explicitly 
and therefore to obtain a locally complete system of exact solutions of the Dirac 
equation with a scalar potential depending on one variable as well as of the 
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stationary Schrodinger equations l|17|) and (|18fl with potentials (|21|l and 12211 

respectively. 

As a first step we construct the adjoint generating pair (see definition [HJ) : 

F* = -/ k and G* = -1. 

Jo 

Next, we write down the expression for the (F, G)-integral: 

Wd {F}G) z = \ (f ( Zl ) Sc / ^\dz - Sc / f (z)W(z)kdz) . 

2 V Jr fo{z) fo{zi) J r J 

By definition, the formal power Z(°)(a, zo; z) for equation H2(J|I has the form 

Z^ia, z ; z) = XF(z)+nG(z), 

where the complex constants A and fi are chosen so that XF(zq) + fiG(zo) = a. 
That is, 

Z (0) (a,z ;z) = \ e P(*)+™*+^y + Me -(P(x)+m a+ iu; 3 /) k 

In order to obtain Z^(a, z ; z) we should take the {F, G)-integral of z[°\a, z ; z), 
where 

z[°\a,z ;z) = \\F x {z) +^iGi(z), 
with XiFi(zq) + fiiGi(zo) = a. Thus, 

Z^(a,z ;z) = [ (AiFi(C) + wGxCO^gjC 



. l| e -P(2;)+m2;+iwy g c /" e -P(x')-mx' -iuy' ^^ e -P(x')-mx' +iujy' + ^ ie P ( x ')+ mx ' ~ iu y' tyd^ 
e -P(x)-ma;-iwj/j i g c / g P(a:')+ma;'+iwi/' j c ^ ie -P(x')-rna;'+iwi,' _|_ ^ e P(x')+mx' -iuy' k)^} 



2 Jz 



where £ = x' + y'k. 

For Z( 2 )(a,2 ; z) by definition E| we have 



Z^(a,z ;z) = 2 Z[ 1} (a, z ; C)d(F,0)C (24) 
where z| 1 ^(a, zo; C) in its turn can be found from the equality 

Zl 1 ' (a, Zq; z) = Z ( 2 0) (a, z : 0^f u G x )C- (25) 

'20 



13 



We note that due to periodicity of the generating sequence containing the gen- 
erating pair 

Zi°\a,z ;0 =zW(a,zo;Q. 

The adjoint pair for G\) necessary for the (iq, Gi)-integral in lt25)l has the 
form 

F* = -e r k and G* = e~ T . 

Thus, 

Zf> (a, z ; z) = 

1 f z , , i , i 

— _r e -P(x)-mx+iwy g c / e P(x )+mx -tujy (^\ e p { x )+mx +iujy _|_ ^-Pfi )-mi -tujy JjW^ 
_ e P{x)+mx-i^jy^ g c /" e -F(x')-mx'+iwj/'] i ^ e P(x')+mx'+iwj/' _|_ ^g^^')-'™'-^?/^,^} 



_ e P(x)+mx-io. ! / kSc A (Ae 2i ^'k-/xe- 2 f p ^') +rox '))rfC}. 

J z 

Substitution of this expression into (|24|l gives us the formal power Z^ 2 \a, zq; z), 
and this algorithmically simple procedure can be continued indefinitely. As a 
result we obtain an infinite system of formal powers which at least locally gives 
us a complete system of solutions of <|2(J[) in the sense that any regular solution 
of H2Ufl can be approximated arbitrarily closely by a finite linear combination 
of formal powers (formula dJl). Moreover, as the corresponding generating 
sequence is periodic, theorem 1141 is valid, and therefore we can guarantee the 
convergence of a Taylor expansion in the formal powers to a corresponding 
solution of l|20[l in some neighborhood of zq. 

A similar procedure works also for equation 119|) . Note that the pair of 
functions Fik = e T k and Gik = — e~ T is a generating pair corresponding to 

As any solution of the Schrodinger equation (|17|l with the potential v\ defined 
by lt2*T]l is the scalar part of some solution of (|2U|) , and any solution of i(T%|) with 
the potential i|22|) is the vector part of some solution of Q20f) . the scalar and the 
vector parts of the constructed system of formal powers give us locally complete 
systems of solutions of <|17|) and l|18fl respectively. 

This last result can also be interpreted in the following way. Consider the 
equation 

-A/ + vf = uj 2 f in Q (26) 

where / is a complex twice continuously differentiable function of two real vari- 
ables x and y, and v is a complex valued function of one real variable x, u> is a 
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complex constant. Suppose we are given a particular solution /o = fo(x) of the 
ordinary differential equation 

-0 + ,/„ = O. (27, 

This implies that we are able to represent v in the form v = p' + p 2 where 
P = /o//o- Then we observe that H26[l is precisely equation (|17|l with to = in 
H21(l . Thus our result means that if we are able to solve the ordinary differential 
equation Q27JI then we can construct explicitly a locally complete system of 
exact solutions to l|2()l) for any ui. For this one should consider the bicomplcx 
Vekua equation l|20[l and follow the procedure described above for constructing 
the corresponding system of formal powers. Then the scalar part of the system 
gives us a locally complete system of exact solutions to 126(1 . 
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